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Abstract 


Define  Yq  to  be  a  geometric  convolution  of  X  if  is  the 

sum  of  Nq  i.i.d.  random  variables  distributed  as  X,  where 
is  geometrically  distributed  and  independent  of  X.  It  is  known 
that  if  X  is  non-negative  with  finite  second  moment  then  as 
p  -*■  0,  Yq/EYq  converges  in  distribution  to  an  exponential  distri¬ 
bution  with  mean  1.  We  derive  an  upper  bound  for  d(Y^),  the 
distance  between  Yq  and  an  exponential  with  mean  Y^,  namely  for 
0  <  p  £  1/2,  d(YQ)  £  cp  where  c  =  EX2/  (EX)2.  This  bound  is 
asymptotically  (p  -*■  0)  tight.  Also  derived  is  a  bound  for  d(Y„+Z) 


where  Z  £  0  is  independent  of  Y 


0' 


2 


1.  Introduction.  If  (X^,  i>_l}  is  an  i.i.d.  sequence  and 

Ir 

is  geometrically  distributed  (Pr(NQ=k)=q  p,  k=0,l,2,...)  and 

Nq 

independent  of  {X^},  then  Yq  =  £  X^  is  called  a  geometric 

1  N 

convolution  of  X.  Closely  related  is  the  random  variable  Y  =  £  X^ 

where  N  =  N^+l,  which  is  also  referred  to  as  a  geometric  convolution. 

Geometric  convolutions  arise  naturally  in  many  applied  probability 
models.  A  recent  paper  of  Gertsbakh  (1984)  discusses  a  rich  variety  of 
applications  in  reliability  and  queues  and  surveys  research  in  the  area, 
most  of  which  was  performed  by  Soviet  authors.  Feller  (1971)  Section 
XI. 6  elegantly  discusses  terminating  renewal  processes,  the  time  until 
termination  being  a  geometric  convolution.  Several  authors  have 
studied  random  sampling  or  "thinning"  of  renewal  processes  which  results 
in  new  renewal  processes  with  geometric  convolution  interarrival  times. 
Jacobs  (1986)  investigates  a  geometric  convolution  in  the  context  of 
combining  random  loads  and  waiting  for  the  stress  to  exceed  a  given 
level.  In  a  G | G | 1  queue  in  equilibrium  the  waiting  time  distribution 
is  a  geometric  convolution  and  has  been  studied  in  this  context  by 
Szekli  (1986)  and  Kbllerstrom  (1976).  Finally,  numerous  applications 
arise  in  regenerative  stochastic  processes.  Consider  a  regenerative 
process  (Smith  (1958))  where  in  each  cycle  an  event  A  may  or  may  not 
occur,  independently  of  other  cycles.  The  waiting  time  for  A  to 
occur  is  then  of  the  form: 

vO 

(1.1)  W  -  l"  Xt  +  z  =  Y0  +  Z 


r.v 
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where  A  occurs  for  the  first  time  in  cycle  N^+l,  is  geome- 

tricially  distributed  with  parameter  equal  to  the  probability  of  A 
occurring  during  a  specified  cycle,  and  Z  is  the  waiting  time  from 
the  beginning  of  cycle  N^+l  until  A  occurs. 

Keilson  (1966)  recognized  the  prevalence  of  (1.1)  and  considered 
the  case  of  small  p.  He  showed  that  if  the  are  non-negative 

with  finite  second  moment  then  W/EW  converges  in  distribution  to 
an  exponential  with  mean  1  as  p  -*■  0.  Thus  the  waiting  time  for 
a  rare  event  (small  p)  to  occur  is  approximately  exponential.  Solovyev, 
(1971)  considered  a  sequence  of  random  variables  of  the  form  (1.1)  in 
which  the  distribution  of  X  varies  with  p  and  obtained  conditions 
for  asymptotic  exponent iality  as  p  -»■  0.  Solovyev  also  obtained  error 
bounds  for  the  exponential  approximation. 

In  this  paper  we  seek  to  bound  the  distance  between  a  geometric 
convolution  of  non-negative  random  variables  and  an  exponential  distri¬ 
bution  with  the  same  mean.  This  problem  is  cited  bv  Gertsbakh  (1984) 
as  being  "of  great  interest  for  engineering  applications".  Defining 
X,  Y  and  Y^  as  above,  q  -  1-p,  u  ■  EX,  ^  “  EX',  ,  -  ..  ,/2.. and 
(t)  =  Pr(Yg  >  O  we  derive: 


(1.2) 


.2.(1+ 2y-l) 

!  q  U  1  F  (t) 

0 


_  E1 

e  ^  ♦  1*  . 

0 


Defining  d(Yg)  as  the  SUP  norm  distance  between  and 
exponential  distribution  with  mean  EY^  *  qu/p,  it  follows  fr 
(1.2)  that  for  0  <  p  _<  1/2: 


om 


The  bound  (1.3)  is  asymptotically  (as  p  -*•  0)  sharp  in  the 
following  sense.  Since  the  problem  is  scale  invariant,  the  true 
upper  bound  is  a  function  of  y  and  p,  say  B(y,p).  We  show 


(Section  4)  that  2yp/(2yp+q)  £B(y,p)  and  thus  B  (y  ,p)  =  2yp  (l+o  (1)  )  . 

Bounds  are  also  obtained  for  d(Yg+Z)  where  Z  is  a  non-negative 

* 

random  variable  independent  of  Yq,  for  d(Y),  and  for  d(Y  )  where 
* 

Y  is  the  stationary  renewal  distribution  corresponding  to  Y  and 


The  current  bounds  offer  improvement  over  those  of  Solovyev 
(1971)  in  that  they  are  derived  under  less  restrictive  conditions, 
require  less  information  about  X  and  Z  to  compute,  and  are  in 
general  tighter.  This  comparison  is  discussed  in  Section  4. 

In  the  above  X  is  assumed  non-negative  with  known  first  two 
finite  moments.  In  Section  3  we  restrict  X  further  and  obtain 
improved  bounds.  Perhaps  the  most  interesting  of  these  results  is 
that  if  X  is  assumed  NBUE  (new  better  than  used  in  expectation, 
defined  in  Section  2)  then  d(Yg)  is  exactly  equal  to  p. 

Our  methodology  is  a  combination  of  reliability  and  renewal 
theory  geared  to  exploit  the  fact  that  Yq  is  NWU  (new  worse  than 
used)  as  pointed  out  by  Daley  (1984)  and  Kollerstrom  (1976).  The 
technique  of  studying  random  variables  through  their  aging  properties 
was  developed  by  several  authors,  most  notably  by  Barlow,  Marshall  and 
Proschan,  and  is  lucidly  presented  in  the  text  of  Barlow  and  Proschan 


(1975). 


2.  Bound  for  General  X.  A  distribution  F  on  [0,”)  is  defined  to 
be  NWU  (new  worse  than  used)  if: 


(2.1)  F (t+x)  >_  F(t)F(x)  for  all  t,  x  >_  0  . 


Similarly  NBU  (new  better  than  used)  is  defined  by  reversing  the 
inequality  in  (2.1).  Thus  F  is  NWU  (NBU)  if  its  survival  distribu¬ 
tion  at  age  0  is  stochastically  smaller  (larger)  than  its  survival 
distribution  at  age  x  for  all  x  >  0. 

Let  {X^,  i_^0}  beani.i.d.  sequence  of  non-negative  random 

variables,  let  N^  be  independent  of  this  sequence  with 

,  n 

Pr(N  =k)=q  p,  k  =  0,1,...  and  define  S  *  T  X.,  and  Y  =  SM  . 

.0  M  ^  n  ~  i  0  Nq 

The  following  simple  but  very  useful  result  is  due  to  Daley  (1984) 

and  Kollerstrom  (1976). 

Lemma  2.1.  Y^  is  NWU . 

Proof .  For  t  >  0  let  Y  be  distributed  as  the  condition  distri¬ 
bution  of  Yq— t  given  Yq  >  t.  Define  N£  =  max{k:S^  <_  t}  and  let 

X^  be  distributed  as  the  conditional  distribution  of  S_,  ,,-t  given 
t  N^+l 

Yq  >  t.  Then: 


Y 

t 


st 

=  Xt  +  Y0  I 


st 


where  X  and  Yq  are  independent.  [j 
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Lemma  2.2,  below,  is  a  slight  generalization  of  a  known  result 
(Barlow  and  Proschan  (1975),  p.  162),  the  generalization  allowing 
for  an  atom  at  zero  which  will  be  required  for  Y^. 

Lemma  2.2.  Assume  that  W  is  NWU  with  an  atom  of  size  p  at  zero, 

oo 

but  with  no  other  atoms.  Let  F  be  the  cclf  of  W  and  M  =  £ 

0 

the  renewal  function.  Then: 


r(t)  >  qe"^"  e-(M(t)-1)  >  e-(M(t)-1)  . 


Proof.  Let  (N^(t)-N^(O) ,t  >0}  be  a  non-homogeneous  Poisson  process 
with  E(N^(t)-N^(0))  =  -Ln(F(t)/q).  This  process  has  its  first  event 
epoch  T^  distributed  as  w|w  >  0.  its  next  interarrival  time 
T^-T^  distributed  as  W-T^jW > T^  and  in  general  its  kth  inter¬ 
arrival  time  T. -T  distributed  as  W-T,  Jw>T,  Since  W  is  NWU 

k  k-1  k-1 1  k-1 

we  have: 


(2.2) 


VTk-i,Tr*‘Tk-i  s-c  w  for  k  -  2 


Next  consider  (N2(t)-N2(0) ,  t  ^  0}  where  N2  is  a  renewal 
process  with  interarrival  time  W.  The  first  event  epoch  of 
N2(t)-N9(0)  occurs  at  ^  W|W  >  0  'v  T^.  Subsequent  interarrival 
times  are  distributed  as  W.  It  follows  from  S^  ^  T^  and  (2.2) 
that : 


(2.3) 


Sk  for  k  _>  1  . 
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Thus  from  (2.3): 

(2.4)  N1(t)-N1(0)\t  N2(t)-N2(0)  . 
Taking  expectations  in  (2.4): 

(2.5)  -Ln(F(t)/q)  £  M(t)-q_1  . 


Finally,  the  result  follows  from  (2.5)  and  the  observation  that 
qe  Pq  £  q(H-pq-1)  =  1.  [] 

p 

Corollary  2.1.  Suppose  that  X  £  0  has  F(0)  <  1  and  p2  =  EX  <  °°. 
Then: 


F  (t)  £  qe 


+-oq_1  (tu_1+2y-l)  -Dq_1  (2y-l)  “t/EY0 

=  qe  '  e 


where  y  =  \ij2\i  . 

Proof .  Case  (i) .  We  first  consider  the  case  in  which  X  has  no  atoms. 

Then  Y„  has  an  atom  of  size  p  at  zero  and  no  other  atoms.  Define 

Gp  as  a  geometrically  distributed  random  variable  with  parameter 

p(Pr(G  =k)=q  p,  k  =  0,l,...).  Define  Ny  and  N  as  renewal 
P  *0  x 

processes  with  interarrival  times  Yq  and  X  respectively,  and 

\  =  enyq’  *X  =  ENx* 

Note  that  N,,  has  1+G  renewals  at  zero,  and  G  renewals  at 
Y0  p  p 

each  renewal  epoch  of  in  (0,=°).  Thus: 


(2.6) 


My  (t)  =  q-1+pq_1  (M^(t)-l)  =  l+pqSl^t)  . 

Alternatively  (2.6)  can  be  easily  proved  using  Laplace  transforms. 
Next,  we  note  Lorden's  (1970)  upper  bound  for  the  renewal  function: 


(2.7) 


Mx(t>  £?+  2y 


The  result  now  follows  from  Corollary  2.1,  (2.6),  (2.7)  and 
-1 


EYn  =  pq  U. 


Case  (ii).  Now  consider  the  general  case.  Define  e^  to  be  uniformly 

distributed  on  (0,e  ),  n  =  1,2,...,  with  lim  e  =0,  e  indeoendent 

n  n  n 

of  X,  and  X  =  X+e  .  Since  X  converges  to  X  in  quadratic  mean, 
n  n  n  ^ 


2  2  bP 

u  =  EX  p  and  y  =  EX  /2v  -*  y.  Define  Y  =  Y  X  , 

n  n  n  n  n  0,n  i  n,i’ 


the 


analogue  of  YQ  with  X  replaced  by  X^.  Then  by  choosing  X^  ..  =  X^.+e 


n,i  in 


we  have  E (Y 


2  r  2 

„  -Y  )  =  E(T  e  .)  ->-0  as  n  Thus  Y_.  converges 

0,n  0  J  n,  l  0,n  6 

in  quadratic  mean  and  thus  in  distribution  to  Y^.  By  Case  (i): 


(2.8) 


F  (t)  »  qe-p<1'1(K1+2V1) 

0,n 


It  follows  by  letting  n  -*■  °°  in  (2.8)  that  the  desired  bound  for 

F  (t)  holds  at  all  continuity  points  of  Y  .  But  since  F  is 
0  U  Y0 
right  continuous  and  the  bound  is  a  continuous  function  of  t,  it 

follows  that  the  inequality  must  hold  for  all  t  >_  0.  Q 

A  non-negative  random  variable  X  with  distribution  F  is 

defined  to  be  NWUE  (new  worse  than  used  in  expectation)  if  u  =  EX  <  00 
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(2.10) 

h*(x)  <  u  l  for  all  x  0 

(2.11) 

st  * 

X  <  X*  . 
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Lemma  2.3.  If  X  is  either  NWUE  or  NBUE  then  P(X*,ye)  <_  Px- 


Proof .  Consider  the  NWUE  case,  the  NBUE  case  being  totally  analogous. 

-1 


First,  defining  A  =  {t:F(t)>_et^  }: 


(2.13)  t ?(X*,ue)  £  sup |G(B)  -  f  y  1e“tU  dtf  = 

n  J  v. 


-1 


-ty'1 

(F(t)-e  Cy  )dt 


From  (2.11),  (2.12)  and  (2.13): 


(2.14)  P(X*,ye)  <  y'1 


-1 


-1 


(G(t)-e~tV  ‘)dt  <  y-1  I  (G(t)-e"tU  )dt  =  p„.  ^ 
A  •'0 


Define  Y  to  be  the  stationary  renewal  distribution  corresponding 


to  YQ  and  Y. 


Corollary  2.2.  If  X  0  with  F(0)  <  1  and  <  °°>  then 


V(Y* ,  (EYQ)e)  <  pq  V 


-1 


Proof.  By  simple  computation  p  =  pq  Y.  The  result  now  follows 


from  Lemmas  2.1  and  2.3.  Q 


Recall  that  d(YQ)  =  P(YQ,  (EYQ)e) ,  the  sup  norm  distance  between 


Yq  and  an  exponential  distribution  with  the  same  mean. 


Theorem  2.1.  If  X  >_  0  with  F(0)  <  1  and  y^  <  °°,  then: 


-1,_  ..  -t/EY 

. . .  -pq  (2y-l)  0  .  — 

(l)  qe  e 


-t/EY 


1  F  (t)  <  e 
0 


°  +  YPq  1 


[2yp,  0  <  p  £  1/2 

f  >  =  i  _! 

[Yq  P,  P  >  1/2 

Proof.  The  bound  on  the  left  of  (i)  is  the  conclusion  of  Corollary 


(ii)  d (Yq)  <_  p  max(2y  >Yq  )  =  ,  _1 

l  Yq”  P,  P  >  1/2 


2.1.  The  upper  bound  follows  from  Lemma  2.1,  (2.11)  and  Corollary  2.2 
by: 
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K: 

I". 


iv 

fry* 


LVj 


#5 


-t/EY  _  -t/EY  *  , 

(2.15)  Fv  (t)-e  U  <  Fv*(t)-e  U  1  V(Y  ,(EY  )e)  £  ypq 


0 


-  Y 


Finally  (ii)  follows  from  (i)  noting  that: 


(2.16) 


-t/EY  _  -t/EYn  -1,»  . 

e  °-F  (t)  <  e  °[l-qe~Pq  (2Y-1)] 


0 


-1 


<  1— q (1-pq  "(2Y-D)  =  2yp  . 


Thus,  d(Y^)  —  max(2YP,Ypq  "S  £  2yp  for  p  £  1/2.  [] 


Corollary  2.3.  Under  the  conditions  of  Theorem  2.1, 


d(Y)  £  pq  ^max(l+Yq  ^,2y~1)  ^  pmax (1+y , 2y~1) . 


-1» 


Proof.  Since  F  (t)  =  q  F  (t)  we  can  multiply  all  three  sides  of 

Y0 


inequality  (i)  of  Theorem  2.1  obtaining  an  upper  and  lower  bound  for 


—  -t/EY 

F^(t) .  Subtracting  e  from  both  F^(t)  and  the  upper  bound 


we  obtain: 


(2.17) 


Fy(t)-e-t/EY  < 


=  pq  1(l+Yq  1)  • 


—  -t/EY  — 

Subtracting  F^(t)  from  e  and  from  the  lower  bound  for  F^(t) 


yields : 
(2.18) 


-t/EY  -  ,  .  ,  -t/EY  "t/EY0, 

e  -F^(t)  £  sup(e  -be  ) 


,  ,  -pq  ^(2y-1) 

where  b  -  e 
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For  b  >_  q,  the  right  side  of  (2.18)  equals  pqP//qe^Y  ^ 
which  is  bounded  above  by  p,  while  for  b  <  q  the  right  side 
of  (2.18)  equals  1-b  which  is  bounded  above  by  pq  ^(2y-1)- 
Thus  d (y)  is  bounded  above  by  the  larger  of  pq  ^(1+yq  ^)  and 
pq  '''(2y-1),  and  the  result  is  thus  proved.  [] 

We  next  seek  bounds  for  dCY^-t-Z)  with  Z  _>  0  and  independent 
of  Yq.  A  few  simple  preliminary  results  are  first  presented. 
First,  for  <  C2>  a  routine  calculus  argument  proves: 


c  /c0-c 

c,  c,  1  2  1 


(2.19)  V(Cle,c2e)  =  (1  -  ^)  (^) 


Next,  note  that  for  any  constant  8: 


<1 - - 

C2 


(2.20) 


o^+e.w^e)  =  t?(w1,w2) 


It  follows  from  (2.18)  that  for  any  random  variable  V  independent 
of  and  W2  that: 

(2.21)  P(WX+V,W2+V)  <_  P(W1,W2)  . 


Thus  F (y)  _>  F[  (t-y)+]-F(t)  for  y,  t  £  0.  Consequently: 

(2.24)  Pr(ce+Z >  t)-Pr(ce >  t)  =  E[F(t-Z)+-F(t) ]  £  EF(Z)  . 

Now: 

(2.25)  EF (Z)  =  E[l-e"C  Z]  =  l-i(c_1)  £  l-e“C  EZ  .  Q 

Finally,  we  need  a  simple  but  useful  result: 

Lemma  2.5.  Suppose  that  X  and  Y  are  both  either  stochastically 
larger  or  stochastically  smaller  than  Z.  Then: 

0(X,Y)  £  max(P(X,Z),P(Y,Z))  . 

Proof .  Immediate. 

The  above  inequalities  now  enable  us  to  derive: 

Theorem  2.2.  Let  X  be  as  in  Theorem  2.1,  and  let  Z  £  0  be 
independent  of  X  with  EZ  <  Define  6^  =  EZ/y.  Then: 

(i)  d (Yq+Z)  £  [2y+<5zq  1Jp,  for  0  <  p  £  1/2. 

(ii)  d (Y  )  £  ypq  1  . 

Proof .  By  the  triangle  inequality: 
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(2.26)  d (Yq+Z )  1  t>(Y0+Z,(EY0)e+Z)4^((EY0)e+Z,(E(Y0+Z))e)  . 

By  (2.20)  and  Theorem  2.1: 

(2.27)  t?(Y0+Z,  (EYQ)e+2)  <_  d(YQ)  <  2YP  >  for  0  <  p  £  1/2. 

By  Lemma  2.5: 

(2.28)  P((EYn)e+Z,(E(Y0  +  Z))e)  £  max[P( (EYQ)e+Z  ,  (EYQ)e) 

P((E(Y0  +  Z))e,(EY0)e)]  . 

By  Lemma  2.4: 

-EZ/EY  -pq_1«z  _L 

(2.29)  P((EY0)e+Z,(EY0)e)  1-e  *l-e  1  pq  5Z  • 

By  (2.19): 

EYn  EZ  -1 

(2.30)  0((E(Yo+z))£>(EYo)e)  <  1___T<1_spq  6Z  . 

Result  (i)  now  follows  from  (2 . 27)- (2. 30) . 

Finally  by  (2.12),  Lemma  2.5,  Corollary  2.2,  and  (2.19): 

(2.31)  d(Y*)  <  max[P(Y*,(EY0)e),O((EY*)e,(EY0)e)] 

-1 

<_  pq  y  • 

Thus  (ii)  holds  and  the  proof  of  Theorem  2.2  is  complete. 
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3.  Improved  Bounds  Under  Additional  Assumptions.  In  this  section  we 
outline  the  improvements  in  the  results  under  various  aging  assumptions 
on  the  distribution  of  X. 

3.1.  NBUE.  Suppose  that  X  is  NBUE  distributed.  Then 
Y  =  Yq  +  X*  <_  Yq  +  X  =  Y,  thus  Y  is  NBUE.  Note  that  Y*  is 
the  stationary  renewal  distribution  corresponding  to  both  YQ  and  Y. 
It  thus  follows  from  (2.10)  that: 

(3.1.1)  -2- <  h(^  <  -B-  for  all  t  >  0  . 

U  —  y  —  qu  — 

Consequently: 

_  £2£  _  _  22L 

(3.1.2)  e  qM  <  F  (t+x)/F(p  <  e  P  for  all  t,  x  >  0  . 

—  Y*  Y*  —  ~ 

Thus  for  p  small,  Y  has  an  approximate  lack  of  memory  in 
that  the  residual  age  distribution  varies  with  t  by  at  most  p 
in  sup  norm. 

Note  that  from  (3.1.1),  (EYg)e  Y*  _<  (EY)e  and  from  (2.11), 

Yc  1  Y*  1  Y-  Since  V(Y,YQ)  p  and  V(  (EYQ)  £,  (EY)  e)  <  p  it 
follows  that: 

(3.1.3)  max(d(Y*,Y0),d(Y*,Y),d(Y*,(EY0)e))  <_  p  . 

Furthermore  by  Lemma  2.3: 

(3.1.4)  d(Y*,(EY)c)  <  Py  *  PPX  . 
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From  (3.1.3)  and  (3.1.4),  using  the  methodology  of  Section  2 
it  is  straightforward  to  derive: 


(3.1.5) 


(3.1.6) 


d(Y.)  <  p 


d(Y)  <  p 
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(3.1.7) 

(3.1.8) 

(3.1.9) 

(3.1.10) 


d(Y  )  <_  pp 


-lx 


d (Yq+Z)  <  p (1+q  6Z) 


Fy(t)  _>  e"t//EY-pPx 


-t/EY 


0  -  ,  .  -1  -t/EY 

£  FY(t)  <  q  e 


Note  that  Pr(YQ=0)  -  p/ (1-qPr (X=0) )  >  p.  It  follows  that  for  any 
non-negative  X  with  finite  mean: 


(3.1.11) 


d(YQ)  >  p  . 


V 

ft 


Thus  (3.1.5)  and  (3.1.11)  show  that  for  F  NBUE: 


(3.1.12) 


d(YQ)  *  p  . 


Finally,  we  mention  that  when  <  p/2  we  can  improve  on  (3.1.6) 

by  using  Daley's  (1986)  bound  for  NBUE  distributions  applied  to  Y. 

This  yields: 


(3.1.13) 


d(Y)  £  /2py  =  /2^ 
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3.2.  NBU.  If  X  is  NBU  then  an  argument  similar  to  the  NBUE  case  of 
Section  3.1  shows  that  Y  is  also  NBU.  Now  YQ  is  NWU,  Y  is  NBU 
and  YQ-t  [Yq  >  t  S-  Y-t  |y  >  t  for  all  t  >_  0.  Thus: 

(3.2.1)  Yq  £  Y-t |Y  >  t  ±  Y 

F  (t+x) 

(3.2.2)  qFy(x)  <_  — — -  <  F„(x)  for  all  x,  t  >  0  . 

?Y(t) 

The  residual  age  distributions  thus  cannot  vary  by  more  than  p 
in  sup  norm. 

Since  Y  is  NBU  we  can  derive  an  analogue  of  Corollary  2.1  for  Y: 

(3.2.3)  Fy(t)  <  ePe“t/EY  . 

Combine  (3.1.9),  (3.1.10)  and  (3.2.3)  to  obtain 


-t/EY 


(3.2.4) 


max(e 


0 


-t/EY 
,e  -pp 


X' 


l  ^(t) 


p  -t/EY 
e  e 


3.3.  NWUE.  Assume  that  X  is  NWUE  distributed.  It  follows  from 
the  argument  of  Section  2.1  that  Y  is  also  NWUE.  Thus  by  Lemma  2.3: 


(3.3.1) 


V(Y  , (EY) c)  £  PPX  . 


st 


Since  Y  <_  Y  ((2.11))  it  follows  from  (3.3.1)  that: 


(3.3.2) 


V.)-.-*/"  I  F  <  pp,  . 


3.6.  DFR.  A  random  variable  on  [0,®)  is  defined  to  be  DFR 
(decreasing  failure  rate)  distributed  if  X-t|X>t  is  stochastically 
decreasing  in  t  0.  Shantikumar  (1986)  recently  proved  that  geo¬ 
metric  convolutions  of  DFR  are  DFR.  Thus  if  X  is  DFR  then  so  are  Yq, 
Y  and  Y*.  Using  the  DFR  property  of  Y  and  Y^  it  follows  from 
Brown  (1983)  p.  422  that: 


(3.6.1) 


(3.6.2) 


*  PPX 

max(d(Y),d(Y*))  < 


-  PPX+1 


d(YJ  < 


0 


-EX_ 


0  —  p  +1  py+q 

0 


A  geometric  convolution  of  DFR  random  variables  arises  naturally 
in  the  study  of  time  to  first  failure  for  repairable  systems  (Brown 
(1984a)  o.  611). 

3.7.  IFR .  Assume  that  X  is  IFR  (increasing  failure  rate).  Then  it 
follows  from  Brown  (1984b)  that: 

2 

(3.7.1)  M ( t )  > 

u 

where  M  is  the  renewal  function  corresponding  to  X.  Then  (3.6.1) 


and  an  analogue  of  Corollary  2.1  yields: 


Since  for  X  NBUE,  0  <_  px  <_  1/2,  we  see  that  (3.7,2)  improves 
upon  (3.2.3). 

3.8.  Pr(X=0)  =  ge (0,1) ,  Known.  If  g  *  Pr(X=0)  is  known,  with 

0  <  S  <  1,  then  an  improvement  in  the  bound  for  d(Y^)  can  be 
achieved.  Define  X'  to  be  distributed  as  the  conditional  distri¬ 
bution  of  X  given  X  >  0.  Then: 

Gp  st 

(3.8.1)  I  X  =  I  X'  where  p*  “  p/(l-gq) 

1  1  1 

(3.8.2)  Yx,  =  (1-B)YX  . 

From  Theorem  (2.1): 

(3.8.3)  d (Yq)  1  2yx,P*  -  (2Yxp)(l-B)/(l-Bq)  . 

The  IMRL  class  is  closed  under  the  transformation  X'  -*■  X  for  all 
0  <  B  <  1.  Thus  (3.5.1),  (3.8.1)  and  (3.8.2)  imply  that  for  X  IMRL : 

(3.8.4)  d(YQ)  <  (yxP> (1-6) / (1-Bq)  . 

4 .  Comments  and  Additions. 

4.1.  Consider  X  ^  Bin(l, a),  i.e.  Pr(X  -1)  -  a,  Pr(X  -0)  -  1-a. 

-  a  a  a 

2  2-1 

Then  •  EX^/2(EX^)  *  (2a)  ,  thus  as  a  ranges  from  1  to  0  y^ 

ranges  from  1/2  to  Thus  all  possible  values  of  y  are  assumed 
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by  the  Bin(l, a)  family.  Now,  let  Yq(cO  be  a  sum  of  ,  X^'s. 

Then: 

(4.1.1)  Pr(Y0(a)-0)  -  ^  -  <^>rp  • 

It  follows  from  (4.1.1)  that  if  B(y,p)  denotes  sup  d(Y^), 
over  all  Y^  with  common  (Y,p)  that: 

(4.1.2)  B(y  ,p)  >_  (q4-2yp) YP  =  2yp(1+o(1))  . 


Thus  as  mentioned  in  the  introduction  the  bound  2yp  for  d(YQ) 
is  asymptotically  (p-*0)  sharp. 

4.2.  Solovyev's  (1971)  bounds  require  existence  and  knowledge  of 

EXm  for  some  me(2,3],  and  of  EZ2,  in  addition  to  our  requirements. 

1/ra-l 

Defining  y  =  [ EX™/ (EX)m]  ,  his  bound  are  0(y  p)  as  P  0  while 

m  m 

ours  are  0(yp)  (Solovyev's  bounds  depend  on  y^/m-2  and  break  down 
for  m  =  2) .  By  the  log  convexity  of  moments  (Marhsall  and  Olkin 
(1979)  p.  74),  Ym  1  Y  for  m  >  2.  The  ratio  ym/y  can  range  from 
1  to  00  depending  on  EX™.  It  is  difficult  to  make  a  comparison  to 
cover  all  possible  cases  but  it  appears  that  the  current  bounds  are 
in  general  tighter.  For  example  Solovyev's  bound  for  d(YQ)  corres¬ 
ponding  to  m  =  3  is  6y-j/y  6  times  as  large  as  ours  as  p  -*•  0. 

4.3.  The  bounds  for  F  (t)  derived  in  Sections  2  and  3  immediately 

Y0 

yield  bounds  for  the  renewal  function  of  a  terminating  renewal  process 
(see  Feller  (1971)  Section  XI. 6). 
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*> 

B 


K: 

S 


p 


4.4.  Defining  A  =  2yp,  our  upper  bound  for  d(Yg)  (Theorem  2.1) 

is  A.  I  conjecture  that  the  best  upper  bound  is  A/A+q.  This  bound 
is  achieved  for  the  Bin(l, a)  family  discussed  in  Section  4.1. 

-G<Pn> 

4.5.  Given  a  sequence  Y.  =)  X  . ,  n  =  1,2,...  where  X  and 

-  n  0,n  L1  n,i  n 

p  varv  with  n,  it  follows  from  Theorem  2.1  that  lim  y  p  *  0  is 
n  '  nrn 

a  sufficient  condition  for  exponential  convergence  of  Y  /EY 

u )  n  u  ,n 

This  condition  is  not  necessary  as  a  slight  modification  of  the  example 

on  p.  874  of  Brown  and  Ge  (1984)  demonstrates.  If  the  X  are  all 

n 

NBUE  then  it  follows  from  p.  872  of  Brown  and  Ge  (1984)  that  a 

necessary  and  sufficient  condition  for  exponential  convergence  of 

Y  /EY  is  lim  p  p  =0  (recall  p  =  1-y  in  the  NBUE  case) ,  and 
n  n  n  n  n  n 

a  necessary  and  sufficient  condition  for  exponential  convergence  of 

Y.  ' /EY  is  lim  p  =  0. 

0 ,  n  0,n  rn 

4.6.  A  simple  argument  is  now  presented  to  show  that  under  very  general 
conditions  geom«  ;;ric  convolutions  are  asymptotically  exponential  as  p  -*■  0. 

Consider  a  random  sequence  {X^  i  _>  1}  which  obeys  the  strong  law 
of  large  numbers  for  ue(0,<3°),  that  is: 


(4.6.1) 


Pr (lim  X  =  u)  =  1 
n 

a-** 


t?: 


i  » 

where  X^  =  —  /,  X..  Define  G  to  a  random  variable  which  is  geome- 

n  r  p 


n  n 


r*  D 

trically  distributed  with  parameter  p.  Consider  Y^Cp)  =  £  Now: 


(4.6.2) 


pYo(p)  -  (PG  )XG  . 

P 


•'  ■  v'» 
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i  * 

ft 


It  follows  from  (3.1.12)  that  d(pG^)  *  p  and  thus  pG^ 

converges  in  distribution  to  an  exponential  with  mean  1.  By 

(4.6.1),  XG  -1+'  u,  thus  pY^(p)  converges  in  distribution 
P 

to  an  exponential  with  mean  y. 

In  the  i.i.d.  case  0  <  EX  <  «  suffices  for  exponential 
convergence  of  Y^.  It  is  not  necessary  that  X  be  non-negative, 

9 

or  that  Gp  be  independent  of  {X^,  i  _>  1),  or  that  EX  <  ®. 

It  is  also  seen  that  a  large  variety  of  dependent  sequences 
lead  to  exponential  convergence  of  geometric  convolutions,  for 
example  stationary  ergodic  sequences  with  0  <  y  <  “.  An 
interesting  problem  is  to  obtain  error  bounds  for  d(Y^)  (also 
d (Y)  and  d(Y^+Z))  for  various  classes  of  dependent  sequences 
{X.,  i  _>  1}. 

If  we  relax  (4.6.1)  to  convergence  in  probability  but  impose 
the  condition  that  G^  be  independent  of  {X^},  then  again 
pYq (p)  is  asymptotically  exponential. 
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